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■^ \ Abstract 

l> 

^T^ ■ It has been revealed recently that, in the scale free range, i.e. from the scale of 

(^ , the onset of nonlinear evolution to the scale of dissipation, the velocity and mass 

density fields of cosmic baryon fluid are extremely well described by the self-similar 
log-Poisson hierarchy. As a consequence of this evolution, the relations among var- 
ious physical quantities of cosmic baryon fluid should be scale invariant, if the 
C^ I physical quantities are measured in cells on scales larger than the dissipation scale, 

regardless the baryon fluid is in virialized dark halo, or in pre-virialized state. We 
examine this property with the relation between the Compton parameter of the ther- 
mal Sunyaev-Zel'dovich effect, y{r), and X-ray luminosity, L:,^{r), where r being the 
scale of regions in which y and Lx are measured. According to the self-similar hier- 
archical scenario of nonlinear evolution, one should expect that 1.) in the y(r)-Lx{r) 
relation, y(r) = 10 ''''[Lx(r)]"(''\ the coefficients A{r) and a{r) are scale-invariant; 
2.) The relation y{r) = 10 '''^[Lx(r)]"'^'' given by cells containing collapsed objects 
is also available for cells without collapsed objects, only if r is larger than the dissi- 
pation scale. These two predictions are well established with a scale decomposition 
analysis of observed data, and a comparison of observed y{r)-Lx{r) relation with hy- 
drodynamic simulation samples. The implication of this result on the characteristic 
scales of non-gravitational heating is also addressed. 
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1 Introduction 



Scaling relation of dimensional quantities is very powerful to reveal the dy- 
namical feature of various physical systems. There has been a considerable 
effort devoting to study the correlations and scaling laws of various observable 
quantities of galaxy clusters. Since virialized self-gravitational system is char- 
acterized by one parameter, mass or virial temperature, one can find a set of 
scaling relations among mass, size. X-ray luminosity, temperature, and Comp- 
ton parameter of Sunyaev-Zel'dovich (SZ) effect if the velocity and mass den- 
sity fields of baryon fluid in cluste rs are assumed to be similar to the virialized 
dark matter halo s (IKaiserl. 1 19861). Observed data of galaxy clusters did yield 



scaling relations (Edge fc Stewartl. Il99ll: iDavid et al.l. Il993 : 



Helsdon fc Ponmad . l2000l : IXue fc Wul . l2000l : ICroston et al. 



WuetaL, 1999 



20051 ). However, 



observed scaling relations generally do not support the predictions given by 
the baryon-dark matter si r nilari ty of virialized dark halos ( iHelsdon fc Ponmanl . 



2000l : iLlovd-Davies et all . I2OOOI ). 



Since Newtonian gravity is scale-free, the self gravitational system of collision- 
less dark matter shows scaling behavior if the power spectrum of initial density 
perturbations is scale-free. These sc aling is regard less of whether the under- 
lying gravitational field is virialized ( Peebles! . 1 19801 ). Thus, if the velocity and 
mass density fields of cosmic baryon matter are given by a similar mapping 
of the fields of dark matter, one may expect the scaling relations of clus- 
ters. However, t he similar mapping assumption is corr ect only in linear regime 
( Bi et al.l . Il992l ). but is baseless in nonlinear regime ( Shandarin fc Zeldovichl . 
19891 ). The nonlinear evolution of cosmic baryon fluid leads to statistically de- 
couple of the fluid from dark matter. The statistical properties of the velocity 
and mass density fields of baryon fluid do show deviation from the unde rlying 
dark matter field (jPando et all . I2OO4J : iHe et all . l200,4 iKim et al.l . I2OO5I ). 



Nevertheless, it has been pointed out by IShandarin fc Zeldovichl (119891 ): the 
dynamics of cosmic baryon fiuid in the expanding universe is scale-free, i.e. no 
preferred special scales can be identified in the range from the onset of nonlin- 
ear evolution down to the length scale of dissipation. It likes fully developed 
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turbulence in inertial range. This idea recently received substantial develop- 
ments. With the hydrodynamic simulation sample of the concordance ACDM 
model, the velocity field of cosmic baryon fiuid i s found to be extreme ly well 
described by She -Leveque' s (SL) scaling formula (IShe fc Levequd . Il994j ) in the 
"inertial range" (JHe et al.l . |2006| ) . The SL formula is considered to be the ba- 
sic statistical features of the scale-free evolution of fully developed turbulence. 
Moreover, the SL formula comes from self-similar log-Poisson hierar chy, which 
is related to the hidden syra metry of the Navier-Stokes equations (JDubruUd . 
I994J : IShe fc Waymird . Il995l ). Very recently, it has been shown that the clus- 
tering of the mass density field of the cosmic bar yon fluid can indee d be well 
described by a log-Poisson hierarchical cascade (JLiu fc Fang . l2008l ). All the 
scaling relations and non-Gaussian features predicted from the log-Poisson 
hierarchy are in very good agreement with the hydrodynamic simulation sam- 
ples. 



These results indicate that, in the scale-free range, the nonlinear evolution 
of cosmic baryon fluid reaches a statistically quasi-steady state similar to a 
fully developed turbulence. For turbulence of incompressible fluid, the fluid 
undergoes a self-similar hierarchical evolution from largest to the smallest 
eddies and flnally dissipates into thermal motion. For cosmic baryon fluid, the 
clustering on different scales can also be described by a self-similar hierarchy, 
and the fluid flnally falls and dissipates into thermal motion. 

This scenario motives us to investigate the scaling relations of clusters from 
the self-similar hierarchy of cosmic baryon fluid. If the observed scaling rela- 
tions come from the self-similar hierarchy, one can expect that 1.) the relations 
of dimensional quantities should be scale-free, i.e. all the scale-dependent co- 
efficients of the scaling relations are scale- invariant; 2.) the relations should 
be held only if the scales of considered regions are larger than Jeans length, 
regardless of whether the underlying gravitational field is virialized, i.e. the 
relations given by cells containing collapsed objects is also available for cells 
without collapsed objects, only if the scale of cells is larger than the dissipation 
scale. 

Other relevant motivation comes from the non-gravitational heating of baryon 
gas of clusters. In order to solve the deviation from the similarity of virialized 
dark halos, various models of n on-gravitational heating and cooling of baryonic 



gas have been p r oposed (e.g. IValageas fc Sim Il999l: iTozzi &: Normanl. 12001 



Voit et al.1 . l2002l : IZhang fc Penl . l2003l : IXue fc Wul . l2003l : iNagai et al.l . 120071 ) 



Since these cooling and heating may introduce characteristic scales, the self- 
similar hierarchy will no longer work on these characteristic scales. Therefore, 
it would be worth to detect the scale on which the above-mentioned two 
predictions to be broken. 



We study these properties with the relation between the Compton parameter 



y of SZ effect and X-ray luminosity L^. The thermal SZ effect is due to the 
inverse Compton scattering of cosmic microwave background (CMB) photons 
by hot electrons of baryo n fluid. The Comptori parameter y depends on the 
press ure of electron gas (IZeldovich fc Sunyaevl . 1 19691 : ISunvaev fc Zeldovich . 



1980') ■ There are many works on the y- Lrr, relation (e.g. Ida Silva et al.l . 12004 



Maughanl . 120071 : iBonamente et al.l . 120071 ). We will, however, focus on the above- 



mentioned two points, which have not yet been addressed right now. 

The outline of this paper is as follows. §2 presents the scaling relations y = 
lO^i"^) L'^i'^^ with observed samples, and shows that A(r) and a{r) are scale- 
invariant. §3 describes the hydrodynamic cosmo logical simulation samples. 
The comparison of the scaling relations of simulation samples with observed 
results is presented in §4. The conclusions and discussion are given in §5. 



2 y-L^ Scaling Relations from Observed Samples 



2.1 Data 



To study the scale free properties, we should find the y{r) - L^^r) relations, 
where y{r) and Ly^{r) are, respectively, the Compton parameter of SZ effect 
and X-ray luminosity measured from regions with s patial scale r. The data 
of X-ray lun iinosity of thes e clust ers are taken from [McCarthy et al.l (120031 ) 
(Xrayl) and iMorandi et al.l ( 20071 ) (Xray2). The X-ray luminosity from area 
on comoving scale r is calculated by 



L,{e)9d9, 



where 6^ = r/[{l+z)dji^{z)] is the angular radius corresponding to the comoving 
scale r, and dA^z) is angular diameter distance. L^{9) is proportional to the 
X-ray surface brightness Sx(0), which can be well fitted by /3-model Sx{9) = 
Sxo[l + {9/9^)^]'-^-^^^^ up to e ~ 10 arcmin. 

Similarly, The mean of y within a region on comoving scale r is given by 



y[rj 



9? 



y{9)9d9. 



(2) 



We wi ll use the SZ effect data from JReese et al.l (120021 ) (SZl) and lBonamente et al, 
(120061 ) (SZ2). The former compiled SZ effects of 18 clusters of galaxies span- 
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Fig. 1. y--Lx relation of observational samples SZl+Xrayl (see Table 1). The Comp- 
ton parameter y are given by average over areas on the comving sizes 0.10, 0.20, 0.39, 
0.78, 1.56, and 3.12 h~^ Mpc, respectively. The solid lines indicate the best-fitting 
for all observational samples. 

ning the redshift range of 0.14 < 2; < 0.78, and the later includes 38 clusters in 
the same redshift range. These data are on angular scales up to ~ 2 arcmin, 
of which the corresponded r is on about the same scale as, or larger than, 
the Jeans length on redshift ~ 0.5. Moreover, the ^-dependencies of y{6) are 
well fitted by /3-model y{9) = yo[l + {e/OcfY^-^^''/^ Therefore, it would be 
reasonable to use the (3 model fitted y{r) to study the y{r) - Ly^{r) relation. 
We will check this point below. 



2.2 Result 



Figure 1 plots the relation oiy{r) vs. L^{r) on scales r=0.1, 0.2, 0.39, 0.78, 1.56 
and 3.12/i~^ Mpc respectively. In this figure the SZ and X-ray data are taken 
from SZl and Xrayl, respectively. The cluster A3 7 is excluded as it shows a 
3-cr discrepancy with the distance- redshift relation (JReese et al.l . |2002| ). Three 
clusters, C10016, A611 and A697, are also excluded due to lacking the data of 
X-ray luminosity, and after all, there are totally 14 clusters used in Figure 1. 



Table 1. r-dependence of a{r) and A{r) 



r{h~^) Mpc 


a 


A 


0.10 


0.32±0.06 


-4.21±0.11 


0.20 


0.36±0.08 


-4.50±0.17 


0.39 


0.44±0.10 


-4.92±0.23 


0.78 


0.51±0.11 


-5.35±0.27 


1.56 


0.56±0.12 


-5.74±0.30 


3.12 


0.61±0.13 


-6.12±0.33 



We make a best-fitting of the y{r)-L^{r) relation with a power law y = 
IqM"^) ]^<^ir) fQj. various scales r as displayed in Figure 1. The coefficients a{r) 
and A{r) are listed in Table 1. Both a{r) and A{r) are significantly depen- 
dent on the scale r. The exponent a{r) increases with scale r, while amplitude 
A{r) decreases with r. If the system is given by a similar mapping of v irialized 
halos, the scaling relation should he y oc LH'^ (ICole k. Kaiserl . llQSSi ). which 
means that a is scale- independent and equal to 0.75. Table 1 shows that the 
values of a on all scales are less than 0.75. Accordingly, the baryon fluid on 
those scales should dynamically deviate from a similar mapping of underlying 
virialized gravitational field of dark matter halos. To check the effect of the an- 
gular scales of ~ 2 arcmin (§2.1), we re-calculate the y{r)-L^{r) relation with 
clusters having dA >1000 Mpc, which with r > 0.5 Mpc for angular scales 2 
arcmin. We find that the coefficients a{r) and A{r) are consistent with Table 
1 within 1-cr range. 



Figure 2 shows a{r) and A{r) as functions of r. Note, both a{r) and A{r) can 
be well fitted, respectively by aa+ba log r and aA+bA log r, and therefore, both 
a{r) and A{r) are scale- invariant. The amplitude lO"^'-'"-' actually is a power 
law of r. This result is consistent with the dynamics of self-similar hierarchy. 
To further test this result, we used other data sets of SZ effect and X-ray 
luminosity. The results are also shown in Figure 2. Although these data sets 
are generally different from each other, all results of a{r) and A{r) can well be 
fitted by the straight line of logr. The fitting parameters are listed in Table 
2. They are the same within 1-a errors. It strongly supports the scenario of 
scale-free dynamics. 
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Fig. 2. The scale-dependence of the coefficients a{V) and A{r^ for various obser- 
vational samples: SZl+Xrayl (solid); SZl+Xray2 (dashed); SZ2-|-Xrayl (dotted); 
and SZ2+Xray2 (dot-dashed). 



Table 2. Fitting Results of a(r) and y4(r). 



sample ^ 



ha 



aA 



SZl+Xrayl 0.52±0.05 0.20±0.07 -5.46±0.13 

SZl+Xray2 0.47±0.06 0.21±0.08 -5.33±0.13 

SZ2+Xrayl 0.49±0.05 0.20±0.07 -5.33±0.11 

SZ2+Xray2 0.50±0.04 0.22±0.06 -5.34±0.10 



-1.27±0.17 
-1.24±0.17 
-1.19±0.15 
-1.22±0.13 



_^Z1 and SZ2 from Reese et al. (2002') and Bonamente et al. (2006') . Xrayl 
and Xrya2 from McCarthy et al. (2003) (Xrayl) and Morandi et al. (2007), 
respectively 

It should also be pointed out that the scales of r < 0.39 h~^ Mpc are actually 
less than the Jeans length of baryon fluid. However, the coefficients a{r) and 
A{r) are still following the self-similar straight lines. In addition, a{r) is seen to 
be less than the value 0.75 from virialized halos. It implies that the dynamics 



of baryon fluid on these scales seems still to be scale free. 



3 SZ effect samples of cosmological hydrodynamical simulations 



3. 1 Simulation 



Before embarking on the numerical calculations, we give a brief summary on 
the baryon fluid when it is in the turbulence-like or self-similar hierarchical 
clustering. The dynamics of growth modes of clustering in an expanding uni- 
verse is sket c hed by a stochastic force driven Burg ers equation (JBerera fc Fang . 
1994 : iJoned . Il999l : iMatarrese fc Mohayaed . l2002l ). The turbulence-like behav- 
ior is due to Burgers' turb ulence, which w i ll be developed when the Burgers 
Reynolds nu r nber i s large (jPolyakovl . Il995l : iBoldyrev et al.l . l2004l : iKraichnanl . 
19681 : iLassigl . |2000| ). A turbulent flow in incompressible fluid consists of vor- 
texes, while the clustering of cosmic matter is irrotational, because the modes 
with vorticity of the perturbed mass density field do not grow. The Burgers' 
turbulence made the initially random field to result in a collection of shocks 
and a smooth variation of the field between the shocks. For cosmic baryon 
f luid, the Burge rs Reynolds number generally is larger in nonlinear regime 
( JHe et al.l . |2004| ). Therefore, the baryon fluid in nonlinear regime consists of a 
collection of shocks with various strengths in both high, moderate and even 
low density areas. The kinetic energy of fluid is dissipated due to the shocks 
on various sc ales. The kineti c energy of fluid is effectively converted into ther- 
mal energy (JHe et al.l . |2006| ). In this context, it is clear that the method of 
hydrodynamical simulations should be capable of capturing discontinuities, 
like shocks, and their effects on energy conversion precisely. 



We will use the Weightly Essentially NonOscillatory (WENO) algorithm, 
which is effective to capture shocks and other discontinuities of the baryon 
fluid, and to give pre cise value of the fluid field between the discontinu- 
ities (JFeng et al.l . |2004| ). This algorithm has been tested with 1.) Shock tube; 
2.) the Sedov- Taylor self-similar blast wave solution, or the Bertschinger's 
similarity solution; 3.) Zeldovich pancake. It is effectively to produce the 
baryo n mass den s ity contour, baryon temper ature contour around massive 
halos (JFeng et al.l . 12004 : iHe et al.l . 12004 . 120051). Some other pr operti es of this 
simula tion algorithm can also be found in iFeng et al.l (12004 ) and iHe et al 
fl2004h . 



The simulations are performed in a cubic box 100^ h~^ Mpc^ with a 1024^ 
grid, and the number of dark matter particles is 512^. The mass of the dark 
matter particle is ~ lO^M©, which corresponds to a density resolution about 
0.01 times of the mean density of intergalactic medium (IGM). The grid size 



is then 10 0/1024 ~ 0.10 h ~^ Mpc. This scale is smaller than the Jeans length 
at 2; < 1 (JBi et al.l . 120031 ) . Therefore, the sample is suitable to describe the 
baryon fluid from the dissipation scale to a few ten h~^ Mpc. We use the 
concordance ACDM cosmology model with parameters fim=0.27, r2b=0.044, 
fiA=0.73, h=0.71, (78=0.84, and spectral index n = 1 and the ratio of spe- 
cific heats is 7 = 5/3 . The transfer function is calculated using CMBFAST 



(jSeljak fc Zaldarriagal Il996l ). 



The atomic proces ses including i onizat ion, cooling and heating are modeled 



as the method in iTheuns et al.l (Il998l ). We take a primordial composition 
of H and He (X=0.7Q, y=0.24) and use an ionizing background model of 
Haardt fc Madaul fl200lh . 



The simulations start at the redshift z = 99, and output the density, velocity 
and temperature fields at redshifts z=2, 1, 0.5, and 0. It is easy for hydro- 
dynamic simulation with Eulerian variable to reach low-density regions. They 
are suitable for a uniform analysis of weakly as well as strongly clustered 
fields of baryon fluid. These samples have been successfully applied to reveal 



the self-s i milar hierarchical behavi or of cosmic baryon fluid (JKim et al.l . 12005 



He et al.l . l2006l : iLiu fc Fang] . 120081). to explain th e transmitted flux of HI and 



Hell Lya absorption of quasars (JLiu et al.l . l2006l ). and to study the relations 



between X-ray luminosity and temperature of groups of galaxies ( IZhang et al. 
2006h . 



As mentioned in §1, the self-similar hierarchical scenario predicts that the 
relation y{r) = 10^*-'"''[Lx(r)]'^'^^) given by areas containing collapsed objects 
should also be available for areas without collapsed objects, only if r is larger 
than the dissipation scale. On the other hand, either the effects of heating by 
injecting hot gas or metal cooling are localized in massive halos. Their char- 
acteristic scales are less than the Jeans length. Therefore, one can expect that 
the y — L relations given by observed samples should be the same as that given 
by simulation samples either with or without considering the localized heating 
and cooling processes. To consider the effect of the me tal abundance, we also 
use a sample with metal cooling of IZhang et al.l (120061 ). For this sample. The 
metal cooling and metal line emission is calculated by the phenomenological 
method: 1) assuming an uniform evolving metallicity Z = 0.3ZQ{t/to), to be- 
in g the present universe age; 2) computing the cooling function using the table 
of ISutherland fc Dopital (119931 ). 



3.2 Samples 



When relativistic corrections is negligible, the Compton parameter y of the 
thermal SZ effect along a line of sight, /, with an angular distance 6 from the 



center of a cluster in the plane of the sky is given by 

yi9) = ^ fn,{l,e)Ul,9)dl, (3) 



where ut is the cross section of the Thomson scattering; ng and Te are, re- 
spectively, the number density and temperature of hot electrons. Since H and 
He atoms are almost fully ionized, we take the electron density Ue = p/^e^ 
where p is the density of baryon gas, pe = 2/(1 + X) with a hydrogen abun- 
dance of X = 0.76. 

Using simulated density and temperature fields, we calculate the parameter 
y{6) with Eq.(3). The mean of parameter y{r) on various scales can be obtained 
using the scaling function of the discrete wavelet transform (DWT) 



where 0j,i(x) is the scaling functions related to cell (j, 1) with comoving size 
100/2-' h~^ Mpc and at position 1 = (l^ , I2, I3). The details of the DWT analysis 



can be found in iFang fc Thewd (119981 ). We take the comoving size 100/2-' h~^ 
Mpc with j =10, 9, 8, 7, 6 and 5, corresponding to scales 0.10, 0.20, 0.39, 0.78, 
1.56 and 3.12 h~^ Mpc, which are the same as that used in the analysis of the 
observed samples in last section. 

The total X-ray luminosity of therm al bremsstrah lung emission from a cell 
(j, 1) is calculated by the same way as IZhang et al.l (12006; ) 



<^J^'' = ''^Tfe;i;^/^''M^'''W^'^- (^' 



where e'^-^(x) is the map of X-ray emissivity, and Vj is the volume of cell (j, 1). 

The DWT decomposition has following advantages. First, the set of scaling 
functions are orthogonal and complete, and therefore, the decomposition of 
Eqs.(3) and (4) does not give rise to false correlation. Second, the DWT cell 
(j, 1) with high mass density can directly be used to identify clumpy structures. 
The DWT-identified cells on scale 1.5 h~^ Mpc is statistically the same as 
clusters id e ntified by traditional method, such as the friend-of-friend algorithm 



( iXu et al.l . Il998l ). Therefore, the DWT variables y^i and {L^)j,\ provides a 
uniform description of the Compton parameter and X-ray luminosity of the 
whole field, regardless of the dynamical details of all cells (j,l). Third, the 



DWT variables can be applied directly to non- Gaussian behavior (JHe et al. 



20061 ; iLiu fc Fang . l2008l ) We will use the Harr wavelet (Daubechies 2) to do the 
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Fig. 3. Compton parameter y vs. baryon density pigm of simulation samples at 
redshift z = 0.5, in which y and pigm are the mean over cells with comoving scales 
r = 0.39, 0.78, 1.56, and 3.12 h~^ Mpc, respectively, pigm is in unit of mean mass 
density pigm of the field. 

calculation below. We also repeat the calculations with wavelet Daubechies 4. 
The statistical features given by Daubechies 4 are basically the same as Haar 
wavelet. 



4 Scaling Relations between the SZ Effect and X-ray Luminosity 



4--1 y-pigin and Lx-pigm relations 



Figure 3 shows the relations between the mean Compton parameter y{r) and 
mean mass density Pigm(^) of cells on scale r for the simulation samples at 
redshift z = 0.5. The comoving scales r are taken to be 0.39, 0.78, 1.56 and 
3.12 h^^ Mpc respectively. Figure 3 has a dark area as a bottom envelop of 
the y-pigm distribution. It gives a tight correlation between y and pigm and can 
be described approximately by a power law oi y cc pl^ for all scales r. This 
relation is basically consistent with the so-called adiabatic 'equation of state' 
T oc p^/^ if considering y oc PigmT. 
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Fig. 4. X-ray luminosity Lx vs. baryon density /Oigm of simulation samples at redshift 
z = 0.5. The Lx is the total X-ray luminosity from cell with comoving scale r = 0.39, 
0.78, 1.56, and 3.12 /i~^ Mpc. pigm are the mean of density within the cell, pigm is 
in unit of mean mass density /Oigm of the whole field. 

In Figure 4, the similar analysis has been performed for the relations between 
the total X-ray luminosity Lxir) and mean mass density Pigm('^)- Clearly, there 
is also a dark area as the bottom envelop of the Lx-Pigm distribution. The tight 
correlation of the dark area yields Lx oc p?„^ for all scales. This is expected if 



considering Lx oc pf T^^"^. 



Therefore, y and Lx given by the tight correlations of bottom envelops in 
Figures 3 and 4 imply y oc L'''^^, which is the same as that given by a virialized 
underlying gravitational field. If we adopt the adiabatic 'equation of state' T oc 
p^/^, the y-Ly^ scaling relation should be y oc L^^'^^ . The observed coefficient 
a are less than 0.75 on all scales. The reason for a < 0.75 is clearly shown 
by Figures 3 and 4. For a given Pigm, the distribution of y is significantly 
scattered from the bottom envelop. For sample at z ~ 0.5, the data points of 
y corresponding to pigm ~ 1 scatter in the range from ~ 10~^^ to 10~^. The 
scattering is due to the heating of Burgers' shocks, which leads to the baryon 
fiuid to be multiphasic, i.e. the relation between temperature and mass density 
can not be described by one polytropic equation. The points of y higher than 
the bottom envelop correspond to state ^ yith temperature higher than tho se 
given by the 'equation of state' T oc p^/^ (JHe et all , l2004l ; iKim et all l2005l ). 
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Fig. 5. The scaling relation between the Compton parameter y{r) and X-ray lu- 
minosity Lx{r). The scales r are 0.78, 1.56, and 3.12 h^^ Mpc. The observed data 
points (squares with errorbars) are from SZl-|-Xrayl and the simulation data points 
(dotted) are taken from sample of redshift z < 0.5. 

The multiphases of baryon fluid are seen for all scales. The scattering in the y- 
Pigm distribution is substantial in high mass density areas (pigm > 1) as well as 
low mass density areas (pigm < !)• The L^-p^g^ distribution (Figure 4) shows 
similar scattering. It can be understood that the Burgers' shock heating is not 
only working in high density regions, but also in low density areas. It leads to 
the deviation of a from 0.75. 



4-2 y-L^ scaling relation 



We now study the scaling relation between y{r) and L^{r). Figure 5 plots 
the distribution of the Compton parameter y{r) against the X-ray luminosity 
Lx(r) for the simulation samples at z < 0.5, in which the comoving scales of 
the cell are taken to be 0.78, 1.56 and 3.12 h^^ Mpc respectively. The observed 
points of t/(r)-Lx(r) are taken from the samples SZl+Xrayl. 

On the scale of 3.12 h^^ Mpc, the simulation data can be fitted by observations 
very well. Using the simulation data with log[Lx(^)/10^^erg/s] > 1, we found 
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that the best fitting scaling relation y = lO'^L^ of simulation sample gives 
a = 0.65 ± 0.03 and A = —6.31 ± 0.05, which are the same as observed 
samples shown in Table 1. For the scale of 1.56 h~^ Mpc, Figure 5 also shows 
the consistency between observed and simulated samples. Since a is less than 
0.75, the scaling relation shown in Figure 5 is dominated by data points located 
above the bottom envelop of Figure 3. In other word, the y-L^ relations are 
mainly determined by the structures which are involved in the evolution of 
the Burgers turbulence. 

On scale 0.78 h~^ Mpc, the simulation sample still shows the same trend as 
observed data, but there are fewer points with log[Lj5f (r)/10'^^erg/s] > 2. This 
phenomenon is more serious on scale 0.39 h~^ Mpc. In this case, the best fitting 
scaling relation yields a = 0.65 ± 0.01 and A = —5.47 ± 0.01, which deviates 
from the observed result a = 0.44 ±0.10 and A = -4.92 ±0.23 (Table 1). The 
deviation on scales < 0.78h~^Mpc is expected. Since the scale 0.78 h^^ Mpc 
is typical of the so-called -R2500 of clusters, and it is less than the virialization 
radius of clusters. Within this scale range, dissipation and non-gravitational 
processes are involved. It is beyond the regime of scale-free evolution. 

Because the observed data contains clusters with redshifts higher than 0.5, 
we also made a comparison between observed data and simulation samples of 
z < 1. The results are displayed in Figure 6, which yields almost the same 
results as Figure 5. We also a nalyzed the simulation samples including metal 



cooling of IZhang et al.l (120061 ). The result is given in Figure 7. It shows the 
effect of metal cooling does not change the feature of Figure 5. Although metal 
cooling may have a big effect on L^ for groups. But it will have the similar 
effect on y. The y-Lx scaling relation still keeps self-similar in the inertia range. 



4-3 Available range of the y-L^ scaling relation 



As emphasized above, the scaling relation between physical quantities of cos- 
mic baryon fluid due to the Burgers turbulence should hold for the entire field. 
Those relations inferred from statistical analysis made in regions containing 
collapsed structures, like clusters, should also be applicable in regions with- 
out that structures. Actually, Figures 5-7 have already shown that the y-L^ 
scaling relation works well for all regions containing strong X-ray emission 
log[Lx(r)/10^3grg/s] > 1. 

Figure 8 presents the y-L^ scaling relation covering a much wider range of 
the X-ray luminosity —3 < log[Lx(^)/10'*^erg/s] < 3, from rich clusters to 
weakly clustered areas. It shows clearly that the scaling relation of Figure 1 
is available in the weakly clustered areas. 



14 



C3) 

° -5 



-6 

-4 

° -5 






o -5 



z<1, r=0.78h"''Mpc 



D„ n 



z<1, r=1.56h"''Mpc 



}^S^* 



t:.B*s* 



z<1, r=3.12h"''Mpc 



e^^ 




1 15 2 2 5 3 

■ log[L^(10^^ergs/s)] 



Fig. 6. The scaling relation between the Compton parameter y{r) and X-ray lu- 
minosity Lx{r)- The scale r are 0.78, 1.56, and 3.12 h^^ Mpc respectively. The 
observed data points ( squares with errorbars) are from SZl-|-Xrayl. The simula- 
tion data points (dots) are taken from sample of redshift z < 1. 

Table 3. Fitting Coefficients a and A for Simulation Data 



range oflog[Lx(^)] 



a 



A 



1 < \og[Lx{r)/l0^^eTg/s] < 3 0.65±0.03 -6.31±0.05 

< log[Lx(r)/104='erg/s] < 3 0.63±0.01 -6.27±0.02 

-1 < log[Lx(r)/10^3erg/s] < 3 0.65±0.01 -6.29±0.01 

-2 < log[Lx(r)/10^3erg/s] < 3 0.65±0.01 -6.30±0.01 

-3 < log[Lx(r)/10^3erg/s] < 3 0.59±0.00 -6.32±0.01 

We did a fitting of y-L^ scaling relation in the simulation sample on the scale 
of r = 3.12 h^^ Mpc and z < 0.5 with different ranges of logLx(^)- The best- 
fitting values of the coefficients a and A for various ranges of logLx(^) are 
listed in Table 3. Obviously, the coefficients a and A are almost independent 
of the range of logLx(^)- The scaling relation, y{r) = 10'^^^^[Lx{r)]°'^^^ is very 
stable within 10^^ < L^ < 10^^ erg s~^. If y and L^ are measured from regions 
on scales larger than dissipation scales, the y-L^ scaling relation still holds 
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Fig. 7. The same as Figure 5, but with simulation sample in which the the metal 
cooling and metal line emission are considered. 



regardless of the dynamical details on the dissipation scales. The wide range 
of Lx is consistent with the scenario of self-similar hierarchical evolution, which 
gives a unified description of the dynamics of clustering on various level in the 
scale- free range. 



Similar analysis was made in simulation samples on scale of r = 1.56h~^ Mpc 
and 0.78 h~^ Mpc at redshifts z = 0, 0.5 and 1 respectively. These scaling 
relations are also valid within the entire range of 10^^ < L^ < 10^^ erg s~^ as 
well. Moreover, from Figure 8, we can see that the scaling relations basically 
are redshift-independent within the range z < 1. Along with the decreasing 
redshift, there are more data points with larger logLx- The y-Lx correlation 
shown in Figure 8 can be seen as a tree with root at left-bottom corner, 
and tip at right-top corner. The formation of clustered objects leads to the 
growth of the tip of the tree along the direction given by the scaling relation 
y{r) = 10^W[L^(r)]"W. 
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Fig. 8. The scaling relation between the Compton parameter y and X-ray lumi- 
nosity Lx of simulation data at redshift z = (top), 0.5 (middle), and 1 (bottom) 
respectively. The scale r are 0.78 (left), 1.56 (middle) and 3.12 (right) h~^ Mpc 
respectively. The solid line is given by the fitting of simulation data. 

5 Discussion and Conclusions 



The scaling studied in this paper comes from the scale-free dynamics of cosmic 
baryon fluid in the expanding universe. That is, there is no preferred special 
scales can be identified in the range from the onset of nonlinear evolution 
down to the characteristic length of dissipation. The clustering evolution is 
described by self-similar hierarchy of baryon fluid. This self-similar hierarchy 
is different from the gravitational self-similarity of virialzed system, which is 
characterized by the mass of the system. On the other hand, the scaling from 
the scale-free dynamics is characterized by the scale range, on which physical 
quantities of the mass and velocity fleld of cosmic baryon fluid are measured. 
That is, the dimensional quantities measured in cells with size larger than 
the dissipation characteristic length should satisfy the same scaling relations, 
regardless whether the cells contains massive collapsed objects. The scaling 
will be broken if the scale is less than the characteristic length of dissipation. 

We demonstrated this scaling with the y{r)-L^{r) relation of baryon fleld. The 
observed y{r)-L^{r) relation can be well reproduced with the hydrodynamical 
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simulation. The important point is that the fitting oi y{r)-L^{r) scahng is not 
based on the identified clusters in simulation samples, but using all cells on 
a given scale r, regardless whether the cell contains rich clusters with strong 
X-ray emission. In other words, the scaling relations y{r) = 10"^^''^[Lx(r)]"'^'''' 
can be used to describe cells containing strong X-ray emission, L^ > 10^^ erg 
s~^, as well as very weak X-ray emission regions L^ ~ 10^° erg s~^. This result 
supports the self-similar hierarchical scenario of the clustering of baryon fiuid. 



The observed y{r)-Lx{r) relations starts to show a deviation from simulation 
results on scales smaller than 0.78 h^^ Mpc, which would be the characteristic 
scale of the dissipation. We calculated the Jeans lengths for each objects used 
in our statistics. All of the Jeans lengths are found to be less than 0.8 h~^ Mpc. 
It is very well consistent with our result. Moreover, this result is also consistent 
with previous studies on the non-gravitational heating of clusters. The non- 
gravitational heating generally is considered to be due to the injection of hot 
gas and energy from SN and AGN. However, the significant effect of injecting 
hot gas and energy by SN is mostly arising from dwarf galaxies, i.e. on scales 
much smaller than clusters. The Lya observations of protoclusters shows that 
the feedback of AGN is not strong enough to heat t he gas of clusters within 



comoving size 0.5 h ^ Mpc ( lAdelberger et al.l . |2003| ). Therefore, the effect of 



the SN and AGN heating would be dramatic only on scales smaller than about 
0.5 h~^ Mpc, but probably cannot heat gas withi n cel l s on scales > 1 h~^ Mp c 



to amount of the order of 1 keV per nucleon (iPen 
Actually we find that scaling relations on scales > 0.78h~^ 
after subtracting the contributions fr om the cent r al par t 



19991 : IWu et all . Il999[ ). 



clusters to both L^ and y, as done in iMarkevitchI (Il998l ): iMaughanl (120071 ). 



Mpc are still held 
< 0.2 h-i Mpc ) of 



As an application, the scaling relations of y{r)-L^{r) would be useful for es- 
timating the contami nation of SZ effe ct on CMB, which is important, espe- 
cially, on small scales (ICao et al.l . |2006| ) . Recent simulation has shown that the 
Planck proje c t wou ld be capable of probing y on the order of y = 10^^ — 10^® 
( Dolag et al.l . l2005l ). It might give a direct test on the universal scaling rela- 
tions of y{r)-Lx{r) given by the self-similar hierarchical evolution. 



Our simulation of the baryon fluid is within the Eulerian framework. One can 
also study the nature of intermittency of fluid with Lagrangian point of view. 
In this approach, the hierarchical clustering can be tracked with Lagrangian 
trajectories. It has been found that the intermittent scaling is related to the 
long time correlations in the particle acceleration, n amely, the random fo rces 



driving the particle motion is long range correlated (JMordant et al.l . |2002| ) . It 



would be interesting to investigate the scaling in the Lagrangian scheme. 
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